Cr Spinor Condensate 



A Biaxial or Uniaxial Spin Nematic 
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We show that the newly discovered ^'^Cr Bose condensate in zero magnetic field can be a spin 
nematic of the following kind: A "maximum" polar state, a "co-linear" polar state, or a biaxial 
nematic ferromagnetic state. We also present the phase diagram with a magnetic field in the 
interaction subspace containing the Chromium condensate. It contains many uniaxial and biaxial 
spin nematic phases, which often but not always break time reversal symmetry, and can exist with 
or without spontaneous magnetization. 



Recently, Tilman Pfau's group at Stuttgart has suc- 
ceeded in condensing a Bose gas of ^^CrQ. ^^Cr atoms 
are spin-3 bosons with electronic spin J = 3 and nuclear 
spin / = 0. In contrast, alkali bosons such as ^^Rb and 
^^Na are spin-1 or -2 bosons with J — 1/2 and / = 3/2. 
Since ^^Cr has a J 6 times larger than the alkali's, its 
dipolar interaction is 36 times larger. Indeed, dipolar 
effects have been observed in Cr 2] , although they have 
not been observed in the alkalis. For this reason, Cr is 
sometimes referred to as a "dipolar condensate. This 
characterization is not quite appropriate. Despite their 
observability, the dipolar energy is no more than a few 
present of the total energy. As we shall see, if the spin 
degrees of freedom of Cr is released, it will have many re- 
markable spin nematic structures stabilized by energies 
much larger than the dipole energy. The Cr condensate, 
in short, is a quantum spin nematic. 

At present, experiments on ^^Cr are performed in mag- 
netic traps which freeze the atomic spins. The spin 
degrees of freedom, however, can be released in opti- 
cal traps. In the case of optically trapped *^Rb and 
^^Na, ground states of different magnetic structures have 
been discovered. Since Cr atoms have a higher spin, the 
number of possible phases will increase. Cr also differs 
from the alkalis in that its scattering lengths in different 
angular momentum channels are very different, whereas 
they are very similar in the alkalis. According to ref.Q, 
the closer these scattering lengths, the weaker the spin- 
dependent interactions. Thus the ratio between spin- 
dependent and density-density interactions in the alkalis 
is about 1% while it is of order 1 in ^^Cr. In other words, 
once the spins of Cr are unfrozen, the spin- dependent in- 
teractions will overwhelm the dipole energy. Any novel 
properties of Cr spinor BEC will be consequences of the 
spin interaction with dipolar effects as a perturbation. 

For spin-3 bosons, the interactions are specified by the 
scattering lengths as in total angular momentum channel 
S = 0,2,4/6. At present, all as except oq have been 
determined[2|- The physical system therefore lies on a 
line in the interaction parameter space, which we refer to 
as the "Cr"-line. While we shall focus on this line, we also 
present the phase diagram as a function of interaction 
and magnetic field, for it will be useful for other spin-3 
Bose systems. As we shall see, the phase diagram is full of 



spin nematic phases which are either unaxial or biaxial. 
These phases typically but not always break time reversal 
symmetry, which can come about with and without a 
spontaneous magnetization Q . The discovery of some of 
these nematic phases will be exciting indeed. 

I. The spin-dependent energy functional: Let 
Cm, m = —3 to 3, be the condensate of a 
spin-3 Bose gas, where n is the density, and C^C = 1- 
As shown in ref.^j, the ground state of a homogenous 
or single mode spinor condensate in zero magnetic field 
is determined by the energy E/Eo = X]s=o 2 4 6 ^-s^S' 
fls = O's/o.B, where Eq = J {2'Kh'^aB/M)n'^, as is the 



Bohr radius, Vs ~ J2m= 



S,M\ 



^ and Bs,M = J2r 



CmCm' (5, M|3, to; 3, m') is the amplitude of forming a 
boson pair in state state \S,M). In the case of spin-1 
alkali gases, (where 5* = 0, 2 in the energy E), the ag's 
differ from each other only by about a few percent, so E 
is almost an identity. In contrast, ^^Cr has a2 = —las, 
04 = bSas, oq — 1 1 2a R, while oq is undermined^]. Us- 
ing the method in ref. a, E can be written as0| 



E = Eo 
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V7B2,M, e = V7Bo,o, 



IBmP and Vn = 7|e|2. The coefii- 



where (S) 

so that V2 = 7J2m \^m\ a^i^^ /-Q 
cients are defined as a = [{qq — Oq) -\- ^^{aQ — 04)], (3 

= [("2 - as) + jf (as - 04)], 7 = jY^iae - 04), and 
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og — Y[{aQ — 04)]. We then have (3 = —4.38, 
7 — 4.91, C = 67.8, and a undetermined. Using the 
relation between 81-82 , (81 • 82)^ and the projection 
operators we can further express eq.lQJ in terms of 



the nematic tensor Af^^ 



E/Eo=-a\Q'^ 



■ ;9Tr7V2 ^ 



V S^S,))/2 as 
7(S>.(S)+C, 



(2) 



= /?/18,7 
-4.38, 7 = 



= 7-^/3, C = C-10/3. 
2.36, C= 111.56 and a 



where a = a — |/3, (3 
For Cr, we have (i = 
is unknown. 

II. Dipole energy : To get an idea of the general 
structure of we use the single mode approximation 
so that all components of 4"^ have the same spatial func- 
tion, ^'^(r) = Y^n(r)C^, and has no spatial depen- 
dence. For harmonic traps Vt = X]j=i "^1^(1" ' ^O^i 



where (ai, a2, 3.3) is an orthonormal triad, the dipole en- 
ergy is (see footnote 7]) 
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i?c/K = ^r7,((S).a,)2, (3) 

1=1 

where ij, = jg^Cg/iB)^ / dridr2n(ri)n(r2)[l - 3(ai • 
r)'^]/r^, r = ri — r2. In particular, we denote the di- 
rection with smaUest 77 as 33. Combining the 7 term in 
eq.(P) with eq.©, all (S)-dependent terms in the energy 
can be written as in eq.© with rji ^ rji + 7. The effect 
of the dipole energy is to align (S) along the direction 
a.3. The phase diagram in zero magnetic field including 
dipole energy can therefore be obtained by first mini- 
mizing eq.QJ (without dipole energy), then replace 7 by 
7 — 773, and align the spontaneous magnetization of those 
phases that acquire it along 3.3. 

To estimate the magnitude of 77^, it is straightforward 
to show that 77, = 4;^?;;°)m ^3^' = where /, = 

/q l"(q)P(l - 3(a, -4)2)//^ |7^(q)p. Approximating n{r) 

™ ^ „^„„„.™. „„ .„^-i^ „ ^ i=i » ^ where Ri is 

the radius of the atom cloud along a^and for cylindrical 
traps {Rx — Ry = R± ^ Rz), Iz varies from to 0.7 as 
Rz/R± varies from 1 to 10. Thus, 77^ is at most about 
25% of 7. As we shall see, all the phases along the Cr-line 
either have zero magnetization, or a weak spontaneous 
magnetization such that the 7 term in eq. ^ contributes 
very little to the total energy. Since Ed is dominated by 
the 7 term, it can be ignored as a first approximation. 

III. Relation between singlet amplitude 8, mag- 
netization (S), and the nematic tensor M: The 
quantities 8, (S), and M provide a characterization of 
the phases. They are, however, not independent quan- 
tities. Let us first consider O, which is 8 = 2C3C-3 ^ 

2C2C-2 + 2C1C-1 - Co' = ^1:1' CCm where A^^^, = 
(— l)"*"*"^ 5m+m',o is the matrix for spin change under 
time reversal. It is also easy to show that the maximum 
value of |8p is 1, and the condition for this is that C 
is invariant (up to a phase x) under time reversal, i.e. 
(A(°)C)m = e^^Cm- This also implies (S) = ^.Thus 
any state with |8| < 1 breaks time reversal symmetry. 

In the presence of a magnetic field B = i?z, the rota- 
tional invariance of eq.l^^) implies that 
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E/E, = a\Q\' + P \Bm?+i{^?-p{Sz) + C, (4) 

A/=-2 

where p — gfisB / Eq. Next, we note that the ten- 
sor Af is hermitian. It then has the diagonal form 
Mi = Aieneij + \2e21e2j + X^e^ie^j, where {6^} are 
the principal axes, (a = 1,2,3), e^ • e^ = 5ab, and the 
= ((Sq • S)^) > are the eigenvalues, satsifying 

Ai + A2 + A3 = 5(5' 4- 1) = 12, 5- = 3. (5) 
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FIG. 1: Phase diagram in zero magnetic field for 7 > 0. The 
dashed line (referred to as the "Cr"-line) is where the physical 
system lies. The solid lines are first order boundaries. The 
nature of the phases displayed are discussed in Fig. 3. 

We shall order the A's so that Ai < A2 < A3. The axes 
(ei, 62, 63) will be referred to as minor, middle, and ma- 
jor axis respectively. Without loss of generality, we can 
always arrange (ei, 62, 63) in a right handed triad. Using 
the terminology of liquid crystals, the system is referred 
as a (spin) nematic if A/" is not isotropic, i.e. not propor- 
tional to the identity matrix. Systems with two identical 
eigenvalues and three unequal eigenvalues will be referred 
to as uniaxial and biaxial nematics respectively. 

The origin of different phases is due the competition 
of 8, (S), and A/", which are not independent quanti- 
ties. Such competition can be illustrated by consider- 
ing /? < 0, in which case TrA/"^ in eq.© favors A3 = 9, 
Ai = A2 = 3/2, which can be achieved by either the po- 
lar state (1, 0, 0, 0, 0, 0, 1) / \/2, or the ferromagnetic state 
(1, 0, 0, 0, 0, 0, 0). However, neither of these states are fa- 
vored simultaneously by the 8 term when a > 0, and the 
(S)^ term, since 7 > for Cr. One of the key features we 
find below is that magnetization (be it spontaneous or 
induced) is generally accompanied with biaxial nematic- 
ity. Even though the energy functional is rotationally 
symmetric about the external field z, this symmetry is 
broken spontaneously in the biaxial nematic state. 

IV. Phase diagram in zero magnetic field, p = 0: 
We have minimized eq.(^ by a combination of analytic 
and numerical methods. The phase diagrams are shown 
in fig.l and 2. Since a is unknown for ^^Cr, the physi- 
cal system lies on the dotted straight line (or "Cr"-line 
for short) in fig.l, — —0.892. It passes through the 
phases A, _B, and C, which we name "maximum polar" 
state (A), co-linear polar state (B), and biaxial nematic 
ferromagnetic state (C), respectively. All order param- 
eters C displayed below are unique up to a phase factor 
and an arbitrary spin rotation. 

Symbols and notations: We characterize each phase 
by its condensate wavefunction C, singlet amplitude 8, 



3 




FIG. 2: Phase diagram in zero magnetic field for 7 < 0. Tfie 
solid and grey lines are first and second order boundaries, 
respectively. The nature of the phases is discussed in Fig. 3. 
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FIG. 3: Phases present at p = 0. a, b, c, d and 5 are real 
numbers. 



magnetization (S) and nematic tensor A/". An isotropic 
N (Ai = A2 = A3) will be denoted by a sphere. Uniax- 
ial nematics with (Ai A2 < A3) and (Ai < A2 = A3) 
are represented by a long and flat cylinder, with major 
and minor axis being the symmetry axis of the cylinder 
respectively. Biaxial nematics with (Ai < A2 < A3) are 
represented as a "brick" with edge lengths given by \a- 
The principal axis Ba is parallel to the edge with length 
\a- Figure 3 shows the structure of the phases. In par- 
ticular, (A) Maximum polar state: This state has time 
reversal symmetry, since = 1. It is a uniaxial nematic 
with A3 = 9, Ai,2 = 3/2. (B) Co-linear polar state: 
It is a superposition of two polar states (1, 0, 0, 0, 0, 0, 1) 
and (0, 0, 1, 0, 1, 0, 0) with unequal weight and a relative 
phase. This system is non-magnetic (S) = and yet has 
broken time reversal symmtry, < |0| < 1. The lat- 
ter is achieved by the phase angle 5, which varies from 
1.3 to 1.4 along the Cr-line from bottom to top. The 
system is uniaxial nematic with Ai < A2 = A3. (C) 



Biaxial nematic ferromagnet : The system has a weak 
spontaneous magnetization (S) = Afz, M = 0.07 to 
0.08, and is a biaxial nematic with either 62 or 63 along 
(S). Since M is small, the contribution of the 7 term in 
eq.([Pl, (and hence dipolar energy) is a very small con- 
tribution to the total energy. This structure is driven 
largely by the competition between the nematic and the 
singlet contribution. The boundary separating B and C 
is a = 3|/3|/(3|/3| +2), 5 = a/7, /3 = P/j. The 
behaviors of Aa and |(S)| are shown in fig. 4. The proper- 
ties of these phases are tabulated in fig. 4. The boundary 
between G and H is a = {2if3+5f3'^)/{36 + 24^+5(3^) 0. 

V. Phase diagram in non-zero magnetic field: 
Figure 5 shows the phase diagram along the Cr-line in 
figure 1 for p ^ obtained by minimizing eq. Q . It has 
an intricate structure near a = 0. All states have non- 
zero magnetization (S) along the direction of the mag- 
netic field z, represented as an arrow. The main features 
of our results are: 

(i) The transition Ai — > Zi — > Gi is a rotation of uniax- 
ial nematic tensor in the Ai phase (represented as a long 
cylinder) along its middle principal axis 62 as one crosses 
the phase Zi , with the rotation angle finally reaching 90° 
at the Gi phase. However, in the Zi phase, J\f acquires 
biaxial nematicity, which continues into the Gi phase. 

(ii) As p ^ 0, for a/7 < -0.372 and a/7 > -0.372, 
the states Cl^ and Cgi reduce to (1, 0, 0, 0, 0, 0, l)/v^ 
and (0, %/3, 0, VW, 0, V3, 0) /4 resp., which are related to 
each other by a 90° rotation about y. (iii) The tran- 
sition Ai —^ Zi —^ Hi is a similar rotation process as 
in (i) except that the rotational angle is zero when one 
reaches the Hi phase. The nematic tensor Af of Hi is 
again uniaxial, except that Hi has zero singlet ampli- 
tude, unlike Ai. (iv) Gi and G2 have similar tensor 
J\f. There is, however, a jump in Q across their bound- 
ary, (v) As the B phase in figure 1 extends into the 
Bi phase in finite field p, its nematic tensor Af changes 
from uniaxial to biaxial, merging into the Ci phase, the 
finite field extension of C. (vi) Among all phases in 



10 





1 






f 






^3! 


^3 






\ - 








1 1 


5 




a/y 





0.08 



0.06 

A 

0.04 V 



0.02 



FIG. 4: The eigenvalues Aa of M and the magnetization (S) 
along the Cr-line in figure 1. 
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FIG. 5; Top panel: Phase diagram for non-zero magnetic field 
along the Cr-line in fig.l. The solid and grey lines are first and 
second order boundaries, respectively. Bottom panel: Expla- 
nation of the phases in the diagram. The numbers a, b, c, d 
are real. All entries of the wavefunctions for the Zj phases 
are generally non-zero. It is not clear they can be represented 
in a simple form as the other phases. 

finite field, Zi,Z2 and Z3 are the ones where none of 
their principal axes are ahgned with z. The Zi has 
it middle axis 62 -L z. In both Z2 and Z3, none of 
the direction cosines ea ■ z vanish. Among these three 
phases, Z3 is the only phase that has vanishing singlet 
amplitude. The phase boundary between Ci — Z3 and 
Z3 — Hi are straight lines at p/7 = 0.667 and 1 resp. 
(vii) The system becomes a full ferromagnet (FF) when 
p/7 > 7.5. To observe the spinor feature, we then need 
B < Be, where B^ = {7.5)-f{2Trh'^aBn/M)/{giiB)- Using 
g = 2,7 = 4.91, we have B^ = [8.34 x 10-^^n]Gauss, 
where n is in units of cm~'^. Thus, with 71 = 5 x 10^^ to 
10"cm-3, Be = 4.2 X 10"^ to lO"* Gauss. 

Detection : Although the ground states of eq.lQ) are 
determined up to an arbitrary rotation about z, this de- 



generacy can be lifted by the anisotropy of a trap through 
dipole interaction eq.Q. With the principal axes (i.e. 
{ea}) fixed by the magnetic field and the trap, the eigen- 
values Xa can be determined by performing Stern-Gerlach 
experiments along the axes Bq. 

Final Remarks: The realization of quantum spin 
nematics, especially the biaxial ones, will be an excit- 
ing development in both cold atom/condensed matter 
physics. These are novel states yet to be realized in 
solid state systems, and biaxial nematics are known to 
have nonabelian defects. Although reducing the mag- 
netic field to the spinor condensate regime for Cr is a 
challenging task, screening a field to 10~^Gauss is not 
outside the reach of current technology. Furthermore, in 
optical lattice setting, one can further increase the den- 
sity with each lattice site, making it easier to reach the 
spinor regime. From the present work, it will not be sur- 
prising if spin nematics can also be found in other higher 
spin Bose gases. 
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